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A NOTE  CONCERNING  A GYROELECTRIC  MEDIUM 


Abstract 


The  fact  that  a homogeneous  electron  gas  when  immersed  in  a 
uniform  magnetostatic  field  becomes  electrically  anisotropic,  i.e., 
gyroelectric,  is  placed  in  evidence.  The  permeability  of  the  gas 
remains  equal  to  .that  of  free  space,  but  its  dielectric  constant  is 
transformed  to  a dyadic  or  tensor  upon  application  of  the  magneto- 
static field.  The  properties  of  the  dielectric  tensor  are  such  that 
a plane  electromagnetic  wave  propagating  through  such  a medium  under- 
goes a Faraday  rotation.  This  rotation  is  the  dual  of  the  Faraday 
rotation  produced  by  gyromagnetic  media. 

The  dielectric  tensor  of  the  electron  gas  is  deduced  and  the 
Faraday  rotation  constant  is  calculated. 

Introduction 

It  is  well  known  that  an  idealized  electron  gas  has  a permea- 
bility u.  and  a dielectric  constant  £ = £-  (I-oj^/o^)  where  u 

and  £q  are  the  permeability  and  dielectric  constant  of  free  space, 
cuQ  is  the  plasma  frequency,  and  cd  is  the  frequency  of  the  wave 
propagating  through  the  gas . ^ ^ 

However,  when  a uniform  magnetostatic  field  is  applied,  the  gas 
becomes  electrically  anisotropic.  This  anisotropy  is  due  to  the  fact 
that  when  an  electromagnetic  wave  propagates  through  such  a medium  the 
convection  current  of  electrons  has  net  the  same  direction  as  the 
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electric  vector  of  the  wave.  Consequently  in  3uch  a medium  B = p H and 

*'"*  o 

D « (£)  • E where  (£)  .is  the  dielectric  tensor. 

The  calculation  of  the  components  of  the  dielectric  tensor  con- 
sists of  first  calculating  the  velocity  of  the  electrons  in  terms  of  the 
electric  vector  E of  the  wave  propagating  through  the  medium,  and  then, 
from  a knowledge  of  this  velocity,  deducing  the  convection  current  of  the 
electrons.  By  adding  to  this  convection  current  the  free  space  displace- 
ment current  of  the  wave,  the  total  current  density  is  obtained.  The  total 
current  density  is  thought  cf  as  a displacement  current  in  an  anisotropic 
dielectric,  and  thus  the  components  of  (c)  are  calculated. 

And  from  the  components  of  ( £ ) expressed  in  terms  of  co,  a>0,  and 

to  where  oj  is  the  gyrofrequency,  in  a straightforward  manner  the 
6 g 

properties  of  a plane  electromagnetic  wave  traveling  through  such  a medium 
in  an  arbitrary  direction  with  respect  to  the  magnetostatic  field  are  de- 
termined. 

The  direction  of  propagation  of  importance  to  microwave  applica- 
tions is  the  one  in  which  the  wave  travels  parallel  to  the  applied 
magnetostatic  field  because  it  is  in  this  direction  that  the  ga3  produces 
a Faraday  rotation  of  the  traveling  wave.  The  Faraday  rotation  factor  is 
easily  determined  from  the  fact  that  the  two  waves  which  travel  in  this 
parallel  direction  are  circularly  polarized  in  opposite  senses  and  have 
unequal  velocities.  And  hence  a superposition  of  the  two  yields  a "linearly" 
polarized  wave  which  rotates  about  the  direction  of  the  applied  magneto- 
static field. 

And,  of  course,  for  propagation  perpendicular  to  the  direction 
of  the  applied  magnetostatic  field  the  Coiton-Mouton  effect^ 


is  exhibited 
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Since  an  electron  gas  with  applied  magnetic  field  produces 

Faraday  rotation,  it  can  be  used  as  a nonreciprocal  microwave  circuit  com- 

(li) 

ponent.  This  aspect  of  the  problem  will  bs  discussed  elsewhere  in 
detail. 


Dielectric  Tensor  of  an  Electron  Gas  with  Applied  Magnetostatic  Field. 

We  deduce  the  dielectric  tensor  of  a medium  consisting  of  a homo- 
geneous electron  gas  immersed  in  a magnetostatic  field  Bq  . When  a plane 
electromagnetic  wave  whose  electric  and  magnetic  vectors  are  respectively 
E and  H travels  through  the  medium,  each  electron  is  subjected  to  a 
force  F which  depends  on  the  electronic  charge  q and  the  electronic 
velocity  v according  to  the  well-known  relation  of  Lorentz: 


F *»  q E + q u vxH  + qvxB 

- - O - “ — -O 


(1) 


wherein  p0  is  the  permeability  of  free  space.  Applying  Newton's  law  to 
each  electron  of  mass  m and  neglecting  the  second  term  on  the  right  side 
of  (1)  we  obtain 


m 


d t 


v ■ qE+qvxB 


(2) 


In  this  equation  v and  E are  real  vector  functions  of  space  and  time, 
and  B is  spatially  uniform  and  independent  of  time. 

Our  first  task  is  to  solve  (2)  for  the  velocity  v . To  do  this 
we  differentiate  (2)  with  respect  to  time, 
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m 


dt 


dc  d 

yi 


(3) 


and  then  postmultiply  vectorially  by  , 


Ij. 


m ~^2  - x -o  ‘ q “^T  -x  £o  £0)x  I. 


»'U) 


n 

Multiplying  (u)  by  q/nT  and  using  a well-known  vector  identity  to  trai  s- 
fo™  the  right  side,  we  get 


2 “4  £ - %:  4z  E x 8 + v • B ) B - SL  (B  • B_)  -Jr  Z • (5) 

m dt^  ° —2  dt  “ —*o  m2  dt  — *~o  — o _2  r 


m 


m4- 


m 


Multiplying  (2)  scalar ly  by  Bq  we  get 


m -f—  v - B “ q E • B , 

dt o ~ -o  ’ 

from  which  it  follows  after  multiplication  by  B q2/m3  that 


” (-Jr*  v • B ) B 

dt  - -o  -o 


(E  • B ) B . 
m3 o -o 


1 cT 

And  operating  on  (2)  with we  get 


m at2 


d3  q &-  K . q d^-  n 

dt3  “ m dt*  - ™ dt^  “ 


(6) 


(7) 


(8) 


We  note  that  the  second  term  on  the  right  side  of  (8)  is  identical  to  the 
term  on  the  left  side  of  ($).  And  the  second  term  on  the  right  side  of 
(5)  is  identical  to  the  left  side  of  (7).  With  these  observations  it  is 
clear  that  (5),  (7),  and  (8)  yield 


-T  I * £ S-f 

dt3  m ” 


q d2^ 

» dt2 


— Ex  B„  + 2^  (E  • B )B 
dt  - - m3 


. (9) 


Since  no  term  in  (9)  contains  a product  of  time -dependent  functions,  we  are 
free  to  restrict  the  time  dependence  to  e ^C01'  by  replacing  B/dt  by  -i«  , 
etc.  ThU3  (9)  becomes 
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\ 2 *5 
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The  term  B )~  is  a scalar  and  equal  to  the  square  of  the  gyrofrequency 
co„j 

(11) 
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( s i ) • < ; 0.~ ) . 
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Without  loss  of  generality  we  introduce,  a rectangular  coordinate 

system  so  oriented  that  B lies  along  the  z-axi.s.  If  we  denote  the  unit 

~ o 

vectors  along  the  coordinates  axis  by  a , a , and  a , we  can  write 

—x  -y  z 

B=aB,E*B=»EB  and  ExB_*aEE-aEB  . And  if 
— o — z o ’ — — o z o — — o — r 


-x  y o —y  x o 


we  define  the  plasma  frequency  os  by 


n 

o 


Nq 

m £ 


(12) 


where  N is  the  number  of  electrons  per  unit  volume,  and  is  the 

dielectric  constant  of  vacuum,  we  can  write  (10)  in  the  following  form: 


Nq  v - -iw 


£ w" 
o o 


CO 


CO 


r 2 

£o  “ojfg 


(j 


ico(co2-co2)  y "y  X 

o 


2 2 
to 

E ) + -?  °— S- 


co2  (co2  - CO  2 ) z zj 
& 

(13) 


Nq  v is  the  convection  current  and  to  it  we  must  add  the  displacement  cur- 
rent -ico  £ E in  order  to  obtain  the  total  current  J,  i.e„, 
o ~ “ 


J » -ico  £ E + Nqv  , ("lU) 

— o ~ — 

The  x,  y,  and  z components  of  J are  easily  obtained  from  (13)  and  (lU) • 
They  are 


oi  c CD 
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o st  ? - \ t — a-  V.  I 


kd  - -?--?)  **  ■ 


[it  ^--r^v] 
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E°1 


It  is  suggestive  to  'write  (15)  in  the  following  manner: 


j - -ioi  £ E - i®  & E 
x xx  x xy  y 


-ioa  £ E - ioi  £ E 
yx  x yy  y 


■ ioi  £ E 
zz  z 


where 


£ = £ (1  ” ~ o — i'  " ^yy 

xx  o 2 ^ 

g 

o ~e  r 

C * r ^ Q & = f 

x>-  -1£o  «<«*-«,*>  > 

& 


£ - £ (l  - rs—) 

zz  o'  2 7 

0) 


fe  can  express  (16)  as  a tensor  equations 


J - -1(0(6)  • E 


wnere 


( £ ) is  the  dielectric  tensor  whose  matrix  is 
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(£) 


^xy 

£ 

yy 
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\ 


\ 

0 

£ I 

zz  j 


(19) 


- g „ and  that 

y* 


It  is  important  to  note  that  £ * £ and  £ 

‘ xx  yy  xy 

the  diagonal  components  are  purely  real  whereas  the  off-diagonal  components 
are  purely  imaginary.  To  place  this  in  evidence  we  let 


xy 


i £ ' 

vy 


and 


yx 


it  • 

yx 


(20) 


where  £ 1 and  £ ' are  purely  real.  Consequently  (19)  becomes 

xy  yx 
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P *Pt 

“XX  * ~xy 


tPA 
\ — ' 


-i£  ' & 

yx  yy 


0 \ 
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(21) 
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where 


£ • 

xy 


2 

° o “g 


co(u^  - ) 


£ • 

yx 


(22) 


When  the  applied  magnetostatic  field  vanishes  (Bq — >0)  we  see  that 

£ • — >G  , £ — >0,  and  the  diagonal  terns  become  equal  to 

xy  * yx 

6.  (1  = w^/oa^).  The  electron  gas  is  isotropic  vhan  B » 0 • 

O o ® 

Wave  Propagation  in  Tensor  Dielectric  Medium 

Now  we  investigate  the  propagation  of  a plane  electromagnetic  wave 
through  a medium  whose  permeability  is  p0  and  whose  dielectric  constant 
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I • 
1 


is  ( c. ) . As  shown  in  Figure  1,  n denotes  a unit  vector  in  the  direction 

of  wave  propagation.  We  let  k 
Jp-  denote  the  vector  propagation  constant , 

and  hence  k ■ n — where  v is  the 
/ " "V 

/ phase  velocity.  Also  wa  let  r be  a 

position  vector,  i.e., 

x r«=ax  + ay+a*.  The  electric 

- -x  -y  ~z 

vector  of  a plane  wave  traveling  in 
uhci  direction  n has  the  form 


Fig.  1.  Arbitrary  direction  n 
of  wave  propagation  in  electron 
gas  with  applied  magnetic  field 


ik  * r -icot 
E e “ — e 


whare  E is  a constant.  To  determine  the  equation  that  E must  satisfy, 
~o  "* 

we  note  that  the  two  curl  equations  of  Maxwell  are 

V x H - -i  o (E)  • E (2li 


VxE-  i«  p H 

— O “• 

where  (£)  is  the  dielectric  tensor  and  the  product  (£)  • E is  the 
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It  follows  from  (2l|)  and  (25)  that 


Vx  S - u‘  n (£)  • E 


Since 


_ _ ik  • r 

V x V x E e “ — 

— n 


ik  • r. 


-k  (k  • e “ ")  + k • k E 


ik  • r 


(2?) 


we  obtain  upon  substituting  (23)  into  (27) 


ik  • r ik  • r 2 . . ik  » r 

k • k E s — -k  ( k * E e “ )Bw  u ( £)  • E e “ 

— u * “ — “O  O — n 


(28) 


Canceling  out  the  exponential  factor,  recalling  that  k « n w/v, 
c2  “ Bu  j and  n • n ■ 1,  we  obtain 


0 0 <=2  e„ 


(2?) 


Without  loss  of  generality  we  rotate  the  coordinate  system  about  the  z-axis 

so  that  9 - tr/2  and  n lies  in  the  yz  plane.  Therefore,  in  (29) 

n - a sin  0 + a cos  0 , n • » E sin  <0  * E cos  0,  and  the 

— -y  -z  ’ — -o  oy  os 

x,y,z  components  of  the  vector  equation  (?9)  are 


(30) 


Eoy  (cos2  0 


,2 


-hZ)  ♦ E 


oz 


(-cos  0 sin  0)  *0 


_ V*"  ^2Z 

0 + E • (-cos  0 sin  0)  + EQa  . (sin2  0 - -£  ) - 0 . 

Since  these  three  simultaneous  equations  are  homogeneous,  for  them  to  yield 
a non-trivial  solution  it  is  necessary  that  the  determinant  A of  the  co- 


efficients vanish: 


V2  *-XX. 
(1  - — — — ) 
cv  '-o 


^ u 


if.  Ja} 

c*  co 


2 v4"  ^ w 

(cos  0 ^ -y  ■ ^ ) (-sin  0 co3  0) 


2 £ i 

ft  vc  ^ZZ  I 

(-sin  0 cos  0)  (sin^  0-  -g  ~g — )| 


In  expanding  this  daterminant  we  find  it  convenient  to  introduce  i;  , 


and  defined 


r _ -I£L  _ ^ 

L1  6n  " * £ » 
o 

£ __  £ 

£2  ‘ T?  * 1 T2  . 


? zz 

3 • T- 
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With  a little  algebraic  manipulation  ( 31 ) leads  to 


tan2  0 * 


fxi  i wv2  i \ 

(■y  “ ~p— ) (~o  - -r~/ 
c2  &1  c2  e2 

/£  iu£  Iri-  + i-h 

c c3  c11  ^ L',l  *2 


This  equato-Jii  is  of  the  second  degree  in  -y  end  always  has  two  real  roots « 

c^ 

When  0=0  propagation  is  along  the  z-axis  and  it  follows  from 
(17),  (32),  and  (33)  that 


vc  1 

— ■ * Ifl"  '■ 

c2  «1 


0)(a>-<ii  ) 
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and 
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..2 


+ i 
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m(.03  ♦ 03  ) 

© 


(35) 


Therefore,  the  two  propagation  constants  for  waves  traveling  parallel  to 
Bq  are  given  by  the  following  two  expressions. 

(36) 


(37) 

Moreover,  when  the  propagation  is  along  the  y-axis , i.e,,  perpendicular  to 
Bq  , 0 is  equal  to  tr/2  and  in  this  case  we  have 

(36) 
(3?) 


The  corresponding  two  propagation  factors  are 
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Up  to  now  we  have  deduced  a formula  (33)  which  allows  us  to 
compute  the  phase  velocity  of  a wave  propagating  in  an  arbitrary  direc- 
tion with  respect  to  the  applied  magnetostatic  field.  We  shall  con- 
tinue by  finding  the  components  of  E and  H and  thus  determine  the 
polarization  properties  of  the  medium. 

Parallel  Propagation 

We  have  seen  that  when  the  direction  of  propagation  is  parallel 
to  the  applied  magnetostatic  field  two  independent  waves  are  possible- j 
one  has  a propagation  constant  k • and  the  other  has  a propagation 
constant  k M , as  given  by  (36)  and  ( 37 ) ♦ The  two  waves  have  no  longi- 
tudinal electric  field,  i.e„,  E - 0,  as  can  be  seen  free,  the  third 
of  equations  (30)  upon  setting  0 = 0.  Nor  do  they  have  a longitudinal 
component  of  the  magnetic  field,  i.e.,  Hz  * 0 , as  can  be  verified  by 
using  (2>)  and  recalling  the  fact  that  Ejj.  and  Ey  are  independent  of 
x and  y . herefore,  the  two  waves  are  purely  THM  . 

One  of  the  waves  is  loft-circularly  polarized  and  the  other  is 
right-circular ly  polarized.  Because  they  have  different  velocities,  a 
superposition  of  the  two  yields  a "linearly"  polarized  wave  whose  polari- 
zation direction  rotates  about  the  z-axis.  This  rotation  is  like  the 
Faraday  rotation  of  optically  active  substances. 


Equations  (30)  for  0=0  reduce  to 


2 € 2 8 

v XX . Vc  XT', 

E (1  - 2 ~£  1 ~ “oy  ^ 2 £ ; 

ox  c o °y  c o 


, v2  ~2  £ 

Eox  ^ 2 • E0y  <1  “ “2  T**0  “ 0 

0 o * C 0 

•u  <*$-£>  • ° 


where  the  two  values  of  v^/e*  are  given  by  (3n)  end  (35)*  The  third  of 
these  equations  shows  us  that  ■ 0 . When  vc/e^  is  given  by  (3U), 

th.9  first  equation  of  (U2)  yields 


v2  ^xy 

7 ^7 


E - ^xx 

oy  x-  ^-gS 

c o 


l*o  fcc 


If  we  had  used  the  second  equation  of  (U2)  we  would  have  obtained  the  same 
result.  The  corresponding  propagation  factor  kQ’  is  given  by  (36).  When 
v2/c2  is  given  by  (35),  we  get  in  a similar  way 


Here  the  propagation  factor  is  k J1  as  given  by  (37)* 

c 

Therefore,  the  electric  components  of  the  two  TEM  waves  propagating  parallel 


to  B are 
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•where  A and  C are  amplitude  constants.  The  corresponding  magnetic 
field  components  are  easily  obtained  by  applying  (2$)  to  (1*5)  and  (1*6) 
thus 
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These  two  waves  (primed  and  double-primed)  whan  cast  into  vector  form  are 
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It  is  clear  that  these  are  two  plane  circularly  polarized  TEM  waves.  Their 
sum  is 
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We  see  that  E • H = 0 and  hence  E and  H are  perpendicular,  as  are 

E* , H'  and  E",  HM  . To  study  the  polarization  of  the  composite  wave  we 

consider  the  ratio  E__/E  . From  (53)  we  know  that 
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If  we  choose  the  amplitudes  A and  C to  be  equal,  we  get 


E 

x 


X . / o 

— - cot  ( 


k.'  - k " 


0 • 


(56) 


Therefore , the  rotation  of  the  resultant  vector  E about  the  z-exis  per 
unit  length  of  travel  is  equal  to  (kQ’  - k11  )/2  . Thin  rotation  is  called 
Faraday  rotation  and  we  denote  it  by  X t i»e», 


(57) 


The  rotation  is  either  clockwise  or  counterclockwise  depending  on  whether 
kQ'  >k0"  or  k^'^k  K . Substituting  (36)  and  (37)  into  (5?)  we  get  the 
Faraday  rotation  as  a function  of  fx*equency: 


(58) 


16. 


Perpendicular  Propagation 

When  the  direction  of  propagation  is  along  any  direction  in  2*0 

plana,  two  independent  waves  are  possible.  One  has  a propagation  constant 

k ' , and  me  other,  k"  , . However,  in  contradistinction,  to  the 

n/2  rr/2 

case  of  parallel  propagation,  only  one  of  the  waves  is  TEM  while  the  other 
is  TM  . 

The  TEW  wave  has  an  electric  field  component  parallel  to  B^_  and 

a.  magnetic  field  component  perpendicular  to  B and  n . The  propagation 

— o ~ 

factor  k*  of  this  wave  is  independent  of  B and  hence  is  the  same  as 

v/2 

that  of  a wave  traveling  through  an  electron  gas  without  any  applied  mag- 
netostatic field. 

The  TM  wave  has  a magnetic  component  parallel  to  . Its 
propagation  factor  k"  ^ does  depend  on  . 

The  TEM  wavs  is  linearly  polarized  parallel  to  Bq  . However,  a 
superposition  of  the  TEW  wave  and  the  1M  yields  a composite  wave  which  is 
elliptically  polarized.  By  this  case  wherein  the  rotation  is  due  to  a 
perpendicular  magnetic  field  the  reader  will  be  reminded  of  the  Cotton- 
houton  effect  of  optically  active  substances. 
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